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OPTIMIZING A LIQUID PROPELLANT ROCKET ENGINE 
WITH AN AUTOM ATED COMB USTOR DESIGN CODE 
— AUTOCOM 

by D. S. Hague, R. H. Reichel, R. T. Jones, and C. R. Glatt 
Aerophysics Research Corporation 

SUMMARY 

A digital computer code, AUTOCOM, has been developed as an aide to 
the liquid rocket engine designer. The code considers the combined 
effects of engine performance, stability, pressure drop, injector 
complexity, chamber length, chamber diameter, and mixture ratio 
characteristics. The code has the ability to automatically define 
the optimal chamber design recognizing these diverse engine charac- 
teristics. An optimum design is generated by means of function 
minimization techniques operating on an engine rating which measures 
the actual engine's payload potential loss from a hypothetical 
ideal combustor which has one hundred per cent of theoretical 
C* performance, infinite dancing rate for all modes of instability, 
zero pressure drop, zero chamber length, chamber diameter equal to 
throat diameter, etc. 

The code is applied to the optimization of an existing engine. 
Payload potential is substantially improved by introduction of a 
series of design perturbations. Computer time required to develop 
the improved engine is four minutes on the CDC 6600 computer. 
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1.0 INTRODUCTION 


In designing a liquid rocket combustion chamber the engineer must 
compromise between characteristics such as performance, stability, 
weight, injector complexity, cost, etc. These engine characteristics 
are not items which are directly controlled by the designer. Instead, 
they are complicated functions of the independent design variables 
available to the designer, for example, injector hole size, chamber 
length, etc. To further complicate the problem, frequently there are 
several techniques that can be used to predict how an engine charac- 
teristic (such as performance) varies with the independent design 
variables . 

If the engine designer had infinite funds and time available to him, 
he could design many combustors with different combinations and per- 
mutations of the various independent design variables. Engine 
characteristics could then be calculated for each design with all 
the available techniques. If the designer had the ability to digest 
all of this information, he could then select the optimum design for 
his particular application. The selection would be made, of necessity, 
on the basis of weighting factors applied to both the various engine 
characteristics and the characteristic values predicted by different 
techniques . 

With limited funds and time, the designer can only examine a few 
designs, and, because he is familiar with only a few techniques for 
calculating the various characteristics, he uses only this limited 
set of techniques to test the acceptability of each design. Using this 
approach some characteristics are never determined until after the 
combustor has been built, tested, and often found unacceptable. For 
example, stability characteristics which are particularly difficult 
to assess frequently result in an unacceptable engine design. Usually 
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the designs selected in a project are those that are very similar to 
designs that have been successful in the past. As a result, a design 
of another group that would be better for a particular application 
is frequently neglected or ignored. Similarly, when trouble is 
encountered during the development phase, changes are made to overcome 
the particular problem using past experience instead of determining 
which variable or set of variables could be used to overcome the problem 
with the least sacrifice to other characteristics. 

The work performed under the present contract was directed to the 
development of a generalized computer program to calculate all the 
characteristics of a given combustor design. The program then uses 
a perturbation technique to determine the changes in the design 
variables that produce the greatest improvement in the rating of the 
combustor design. The program then follows the path that produces 
the greatest improvement in the rating to arrive at a combustor design 
that has the best combination of all variables . This design is called 
the optimum combustor design. The automated combustor design code 
which generates the optimum combustor design has been given the acronym 
AUTOCOM. 

In any optimiz .tion situation, the engineer/designer is ultimately faced 
with the problem of selecting the rating or value function which is to 
be minimized. In this report the rating of a design is based on a 
weighted average of all the characteristics of a given combustor. The 
weighting factors are constants used to obtain both average character- 
istics and a rating. The constants are intended to allow the designer 
to introduce his views regarding the importance or validity of one 
technique for obtaining a given characteristic versus another technique 
for obtaining the same characteristic. For example, if the designer 
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believes that only one technique is valid for predicting the perfor- 
mance of a given design , he vill assign a unity value to the weighting 
factor constants for that specific characteristic, and the constants 
for all the other performance characteristics will be zero. Similarly, 
the weighting factor constants in the equation to obtain a single 
rating for a given combustor design are intended to allow the designer 
to introduce the relative importance of different types of character- 
istics, for example, stability versus performance. The weighting factor 
constants, therefore, give the designer the same control and flexibility 
in the computer program as he has in the pi^esent tr cut-and-try n system . 

To establish a base point for the rating system, a hypothetical ideal 
combustor is given a rating of zero. The hypothetical ideal combustor 
would have one hundred per cent of theoretical C* performance, infinite 
damping rate for all modes of instability, zero pressure drop, zero 
chamber length, chamber diameter equal to the throat diameter, etc. 

Specific techniques for obtaining the various combustor characteristics 
contained in the AUTOCOM code are outlined in Appendix A of this report. 
The code is written in a modular fashion which permits rapid extension 
of the combustor characteristic equations. This approach leads to an 
open ended code capable of future development and extension consistant 
with the growth of capability in combustor design analysis. 

The optimum combustor design procedure is now an operational tool capable 
of rapid application to practical design problems. This report is 
primarily intended as a demonstration of the current version of the 
AUTOCOM code. An existing liquid propellant rocket engine having a well 
established rating value is studied. An improved design is then auto- 
matically generated by the AUTOCOM code, and a significantly better design 
is developed. The approach followed is outlined in Section 2; Section 3 
describes the nominal engine in detail. Section H traces the development 
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of the improved design. Conclusions are presented in Section 5> and 
a self-contained brief outline of the AUTOCOM analysis procedure is 
presented in Appendix A. Appendix B describes a recently developed 
multivariable optimization algorithm which is believed to represent a 
significant improvement over other existing algorithms in terms of 
the number of design perturbations required to obtain an optimal design. 
Appendix C presents a list of weighting factor constants used in the 
development of a combustor rating or value function for the study of 
the sample engine; the application of these weighting factor constants 
is discussed in Section 2 . 


2.0 APPROACH 


The AUTOCOM code considers the following major characteristics 
in the combustor design synthesis: 


I Performance 

II Stability 

III Pressure Drop 

IV Injector Complexity 

V Chamber Length 

VI Chamber Diameter 

VII Mixture Ratio 


In view of the uncertainties associated with prediction of combustor 
design characteristics, each major characteristic is computed as an 
average engine characteristic . Each average characteristic is an 
appropriate weighted sum of the characteristic value obtained by altern- 
ative accepted computation procedures. Each such computation procedure 
defines a specific combustor characteristic . The weighting factors 
employed in combining a subset of the specific engine characteristics 
into a particular average engine characteristic may be selected by the 
user. They thus can be used to reflect user relative confidence in 
each specific combustor characteristic. 

Each specific engine characteristic is a function of the design variables 
entering into the combustor design procedure. These combustor design 
variables include 

1. Diameter of fuel orifices 

2. Diameter of oxidizer orifices 

3. Number of fuel orifices 

U, Number of oxidizer orifices 


6 


5. Volume of fuel manifold 

6. Volume of oxidizer manifold 

7. Length of fuel orifices 

8. Length of oxidizer orifices 

9. Length of chamber 

10. Diameter of chamber 

11. Mixture ratio 


A subset of these variables defines each specific engine characteristic. 


The combustor rating provides a single numerical measure of the com- 
bustor’s capability and is constructed on the basis of a weighted sum 
of the average engine characteristics. The weighting factors employed 
in computing the combustor rating are user-defined in the AUT0C0M code. 
In this note, these weighting factors are based on the impact of each 
average engine characteristic on vehicle payload capability; they define 
the payload penalty associated with each characteristic. 


The rating function employed in the AUT0C0M code is 
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where 

F is the average performance characteristic "based on C* 

efficiency and varies from 0 to 100 per cent. 

is the average stability characteristic based on an 
equivalent damping rate and varies from -«> (damps at an 
infinite rate with time) to + «>( grows at an infinite 
rate with time) 


III 


is the average pressure drop characteristic based on the 
pressure drop across the injector face and varies from 
0 tO oo. 



is the average injector complexity characteristic based 
on the number of injector elements, type of element, injector 
cavity volume and injector face thickness: varies from 

0 to 


~F is the average length characteristic based on the chamber 

length from injector to nozzle throat, varies from 0 to oo. 

Fyj is the average chamber diameter characteristic based on 

the chamber diameter at the injector face, varies from 
0 to co . 

average propellant mixture ratio characteristic which 
varies from 0 to ». 


' ’ B rai ’ c fii 


and the constants A pI> B pI ,B FII ,B piII . . 

are weighting factors used to define appropriate measures for combining 
the average engine characteristics into the final combustor rating. 
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The average engine characteristics in turn are appropriate weighted 
averages of the specific combustor characteristics which are computed 
from well-defined equations and/or curves accepted by the engineering 
and scientific community. Weighted averages employed in the AUTOCOM 
code are 
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Here, the specific combustor characteristics, f p , are obtained as 
follows , with definitions given in Appendix A. 
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1*11 is the percent mass vaporized of fuel. 

ff2 is the percent mass vaporized of oxidizer. 

f^^ is the C* efficiency determined by the mixing model of NASA. 

f ^ is the C” efficiency determined by the A. D. Little Correlation 
for Pulsed Combustors. 

f^t- is the C * efficiency determined by the A. D. Little Correlation 
for Non-Pulsed Combustors. 

f U is the chugging decay rate based on the fuel system. 

f 1 is the chugging decay rate based on the oxidizer system. 

b 

f 22 is the stability characteristic based on the A. D. Litule 
Correlation for Pulsed Operation. 

f^^ is the stability characteristic based on the A. D. Little 
Correlation for Won-Pulsed Operation. 

f ^ is the stability decay rate characteristic based on the 
stability analysis of Dykema for the fuel. 

f is the stability decay rate characteristic based on the 
stability analysis of Dykeina for the oxidizer. 

fgg is the stability decay rate characteristic based on the 
sensitive time lag model for a longitudinal mode. 
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is the stability decay rate characteristic based on the 
sensitive time lag model for transverse modes. 

is the stability decay rate characteristic based on the 
response function approach of NASA Levis Research Center. 

is the stability characteristic based on the non-linear 
stability analysis of MSA Levis Research Center. 

is the injector fuel pressure drop characteristic. 

is the injector oxidizer pressure drop characteristic. 

is the number of injector fuel plus oxidizer holes characteristic. 

t 

is the volume of the injector oxidizer dome characteristic, 
is the volume of the injector fuel dome characteristic, 
is the length of the injector oxidizer holes characteristic, 
is the length of the injector fuel holes characteristic, 
is the injector type complexity characteristic^ 
is the injector length characteristic. 

is the combustion chamber diameter characteristic. 


is the propellant mixture ratio characteristic. 



The combustor design optimization process is based on minimization of 
the combustor rating and, hence, the payload penalty. The rating is 
clearly a function of the combustor design variables, and the weighting 
factors entering into both the rating equation and the average engine 
characteristics. In a given computation, these weighting factors are 
fixed, based on payload impact and degree of confidence in each 
specific combustor characteristic. It follows that the combustor 
optimization problem can be formally stated as 

<j>* = Min )] (9) 

where <f> is the combustor rating, <p * is the optimal combustor rating, 
the ctj_ are the combustor design variables, and a± is the the vector of 
these design variables. Equation ( 9 ) defines a multivariable optimi- 
zation problem which, due to the non-analytic nature of several specific 
combustor characteristics, can only be solved by numerical methods. 
Reference 1. These methods involve repetitive combustor design eval- 
uations using perturbed sets of combustor design variables. By properly 
organizing the design variable perturbations on the basis of their 
effect on the combustor rating, the succession of designs generated can 
be made to converge to the optimal design satisfying Equation (9)« 

Selection of successive design variable perturbations involves the appli- 
cation of multivariable search techniques. A variety of such search 
techniques have evolved in recent years. They include elemental one- 
parameter-at-a-time methods, organized methods which require the eval- 
uation of first- and second-order partial derivatives 8<|>/3oti and 
3 2 <j)/8a^8a. and finally randomized techniques. The AUTOCOM code contains 
a selection of all three types of search procedures based on the Ref- 
erences 1 and 2 optimization program AESOP. The searches may be used 
separately or in combination at the user's option. Usually a combination 
of searches will provide more rapid and regular convergence to the optimal 
design than will the repetitive application of a single search algorithm 
such as, for example, steepest-descent . 
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An overall schematic diagram of the AUTOCOM program is presented in 
Figure 1. The remainder of this note describes the application of 
AUTOCOM to the optimization of a liquid rochet engine combustor. The 
procedures employed to insure an adequate numerical model of the 
design process vhile controlling elapsed computer time are described 
in some detail. Convergence from an initial nominal design to the 
final optimal design is reported and convergence plots for the combustor 
rating and each design variable are supplied. An outline of the 
available specific combustor characteristic computations is presented 
in Appendix A. 
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FIGURE 1. OVERALL SCHEMATIC OF AUTOCOM PROGRAM 
































3.0 THE NOMINAL ENGINE 


3.1 DEFCHIPTION 

T'ne nominal engine is an existing 15000 lb f- thrust liquid rocket engine. 
Combustor design variables for this engine are 


1. 

Fuel Orifice Diameter 

.195 inches 0.D.,.1^5 

inches I.B. 



(.129 inches diameter 

equivalent hole) 

2. 

OX Orifice AP Diameter 

.084 inches 


3. 

Number Fuel Elements 

23.6 


4. 

Number OX Elements 

216 


5. 

Vo 1 . Fue 1 Man i f ol d 

3 

10 inches 


6. 

Vol . OX Manifold 

3 

10 inches 


7. 

Length Fuel Orifices 

.06 inches 


8. 

Length OX Orifices 

.40 inches 


9. 

Length of Chamber 

11.22 inches 


10. 

Diameter Chamber 

10.28 inches 


11. 

Mixture Ratio 

5.06 lbm OX/ lbm Fuel 



Other pertinent but fixed design 

Propellants 
Element Tj r pe 
OX Orifice Velocity Diain. 
Element Impingement Angle 
Total Propellant Flow 
Thrust 

LOX Temperature 
H 2 Temperature 
Hi ro at Di amet e r 


parameters include 
Hydrogen and LOX 

Concentric Tube (Hydrogen on outside) 
.11 inches 
0 de gre e s 
33. 72 lbm/second 
15000 lbf . 

Boiling 
3*+9°R 

5.14 inches 
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Tn ■*. s en r t i ne runs 


the foilcvinr measured conditions: 


Fuel Injection AP 
LOX Injection AP 
C* Efficiency 
I 

Sp 

Chamber Pressure (injector 
face / static ) 

These running conditions were used to 
in the combustor synthesis. 


83.1 psi 
)4 8.3 psi 
98.6 per cent 
444 Ihf-sec/lbm 

39o. I; psi 

check the nominal engine description 


Computed running conditions from the AUTOCOM code were 


Fuel Injection A? 

LOX Injection AP 
C* Efficiency 

-^sp 

Chamber Pressure (inj 


82.9 P si 
48.5 P si 
See Section 3.2 
444.8 Ibf-sec/rbm 


static) 


372.8 psi 


It is assumed that the low chamoer pressure c onput e q results from the 
ratio of specific heats employed for the propellant combination (y = 1.2505) 
and the combustion temperature (T c = 5722°R). A more complete set of 
computed engine running conditions is presented in Table I. 


3.2 SPECIFIC CC-: 3' 
Specific combustor 
in Table II. The 
selected specific 
rating value resul 
we i ghr. ;i g factors 


FETOR CHARACTER! STICS 

characteristics for the nominal engine are presented 
average engine characteristics resulting from the 
combustor characteristics are shown in Table III. The 
ring from the selected average engine characteristic 
is also presented in Table III. 


It should be no 
con s u; . 0 s a l.mc s t 


.cd tr.r.« the specific combustor stability characteristic 
ell the computational time required for the evaluation of 
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Tit) tuG]v( TYPE NO.l 

Ir.ST C" Sf 

Tn^uST = 1SUOO POUND 
O.Y / HP PROPELLANT 
OPT J M I 7 A T ION PUN 

GENERAL ENGINE PARAMETERS 
COM-ilSTOP THRUST FORCE 
P d OPFLLAMT SPECIFIC P^ULSE 

chav-rep pressure. at injector head 

TCTAt PROPELLANT FLOW PAT £ 

FUEL WEIGHT FLOW »ATE 

OXIDIZE- 1 WEIGHT FLOW PATE 

COMBUSTION TFMPEPATUWF IN CHAMBER 

PATIO OF SPECIFIC HEAT OF COMBUSTION GAS 

GAS CONSTANT OF COMBUSTION GAS ( ft . lb . /Tb . °R ) 

IDEAL THPUST COEEFICIENT 

ACOUSTICAL LENGTH OF CHAMBER 

MOLECUL'P WEIGHT OF COMBUSTION GAS 

MEAN PE si DEUCE TIME OF GAS IN CHAMBER 

SPEED OE SOUND IN CHAMBEP 

COMPOST I ON CHAMBER MACH NUMBER 

INJECTOR PRESSURE DROP E OR EOEl 

I NJECT0 0 PRESSURE DROP EOR OXIDIZER 

COMBUST I ON ChAmPER VOLUME 

AVERAGE VELOCITY OF GASES IN CHAMBER 

FUEL INJECTION VELOCITY 

OXiniZEP INJECTION VELOCITY 

TOTAL A«FA OF FUEL INJECTOR ORIFICES 

TOTAL AREA OF OXIDIZER INJECTOR ORIFICES 

Chamber LENGTH TO V A P 0 R 1 7 1 50 PER-CENT OE FUEL 

CHAMBER LENGTH TO VAPORIZE 50 °£P-CENT OE OXIDIZER 


THWUST = 

15000.0 

ISP 

A A A • 6 A 0 

PCI = 

37?. 765 

WT 

33.7200 

w r = 

5.S6AA5 

wox = 

26.1556 

TCOMH = 

576A.91 

SPHEAT = 

1 .2A9AA 

R6AS a 

127.127 

CFIDEL = 

1.93928 

LOOM = 

8.8157A 

MOLWT = 

12.1A53 

THETAG = 

8.567267F 

CS 

65133.9 

MC = 

.201068 

OELPF = 

82.90A9 

DELPX = 

A8.A525 

VC 

681 .565 

VELC a 

13096. A 

VFUEL = 

13013.0 

VOX I D = 

339.511 

TAF a 

2.831 BA 

TAA = 

1.19702 

L5FUFL*= 

0. 

L50XI0 = 

.636932 


* Fuel in a gaseous state 

TABLE I. NOMINAL ENGINE RUNNING CONDITIONS 
(Dimensions ore in inches, lb., sec., °r) 
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NUMPFR 3)6 


ENGINE TYPE uO.l 
It ST fA>t 

m-just = 1 1 ><» d o pound 

0 ? / HP PROPELLANT 
GPTIM/4TI0M PUN 

SPfCIfiC COMBUSTOR CHARACTERISTICS 


PER CEJT HtsS FUEL VAPOR IVED Fll 
PFR CENT MASS OF OXIDIZE* VAPORIZED Cl? 
C° EFFICIENCY MIXING MODEL Cl 3 
C» FCFICIENCV Piit.SFD COMBUST OPS CIA 
C° F CP I C I EVCY M 0 'T-PULSFD COMhUSTORS CIS 
FUEL SYSTEM CHUGGING DECAY RATE C 20 
0 ZIDI 7 PR System CHUGGING DECAY RATF. C 21 
PULSED INSTABILITY CHARACTERISTIC C 22 
NON-RULSED INSTABILITY ChaRACTE->IST IC C 23 
DYKEMA FUEL STAFRITY OECaY KATE C 2 A 
DYKFMA OXIDIZER STABILITY DECAY RATE F 25 
STABILITY LONGITUDINAL TI"E LAG F 26 
STABILITY TRANSVERSE TIKE LAG C 27 
STABILITY LUC orSPONSE FUNCTION C 28 
STABILITY PR I EM LINEAR ANALYSIS C 29 
FUEL PRESSURE DROP CHARACTERISTIC C 31 
0 XIDI 7 EP PRESSURE DROP CHARACTERISTIC C 32 
FUEL PLUS OXIDIZER HOLES CHAKATERIST 1 C CA 1 
0 XID 17 EP DOME VOLUME CHARACTERISTIC CA 2 
FUEL DOME VOLUME CHARACTERISTIC CA 3 
Oxini 7 FP HOLE LENGTH Char AC TER I ST I C F 4 A 
FUEL HOLE LENGTH CHARACTERISTIC ^65 
INJECTOR TYPE COMPLEXITY CHARACTER 1 ST IC F 46 
INJECTOR LENGTH CHARACTERISTIC F 51 
CHAMBER DIAMETER CHARACTERISTIC E 61 
MIXTURE PATIO CHARACTERISTIC F 71 


TABLE II. NOMINAL ENGINE SPECIFIC COMBUSTOR CHARAC' 


= 100. (Gaseous) 

= 99.8122 
= 100.000 
= 91.3821 
= 73.7665 
= - 1036 . 

= -26l . 5 

= -37 .4819 

= Not Computed 
= Not Computed 

= - 2964.82 
= -1681 .31 
= -412.643 
= -360.322 
= -46427.6 
= 82.9049 
= 48.4525 
= 431.000 
= 10.0000 
= 10.0000 
= .400000 
= 6 • 0 OOOOOE-02 

= Not Computed 
= 2.18268 
= 2.00000 
= 5.05990 

'ERISTICS 
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a combustor. This is due to the time consuming complex characteristic 
equation solutions required for chugging (f 20 and f 2 l)> longitudinal 
time lag stability analysis (f 26 )> transverse time lag stability analysis 
(fpjK and the Levis response function stability analysis (fpg)* Table 
IV presents a summary of the characteristic stability equation roots for 
the nominal engine. The least stable root is obtained from the trans- 
verse time lag analysis (fp 7 ) using a value of S v ^ = 3 . 05 ^ 3 . 

3.3 A MOTE ON STABILITY ROOTS 

Some difficulty was initially experienced in computing the nominal engine 
stability characteristics for the time lag analyses. The AUT0C0M 
program assumes a value of the Reardon interaction index* n* of 0.5 
for the? longitudinal time lag analyses (f 2 g) an d 1-0 for the transverse 
time lag analysis ^ 27 ). With these interaction index values* the nomi- 
nal engine was found to be slightly unstable in tvo of the transverse 
modes * Table V. A sensitivity study on the effect of interaction index 
value v as undertaken; as a result* an interaction index value of 
0.45 vas subsequently utilized in all transverse time lag analyses and 
an ijiut?racuiun ini .a value of 0.9 vas used for the longitudinal analysis. 
These values were used to obtain the time lag analysis roots shown in 
Table IV. 

A second point should be noted regarding the stability roots. No root 
is found corresponding to the third value of S v ^ = 3.8317 in the transverse 
time lag analysis ^ 27 ). The missing root can be found by varying the 
initial guess value in the complex plane for this particular root. 

Following this procedure the missing root vas found to be at the point 
(-.59^2^ - J3.51W) when the Reardon interaction index vas 0.9* The 
root is thus highly damped. 
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When coi; fencing an optimization study, side analyses of the above type 
may be required to locate particularly difficult roots. This procedure 
should also be followed whenever the root imaginary part is not 
approximately equal to the corresponding value of Sv^ in the transverse 
time lag analysis (f 27 ) and the Bessel argument, m, in the Lewis 
response function analysis (f28)- This point is discussed further in 
Section U.l. 
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1* .0 OPTIMIZATION COMPUTATIONS 


Optimisation computations were initially undertaken using all specific 
combustor characteristics and all stability roots. However, it was 
noted that the design variable perturbations introduced little change 
in the computer time consuming stability equation roots. Accordingly 
the combustor analysis was divided into two classes of computation. 
These were an approximate analysis which considered fewer (possibly 
none) of the stability roots and a complete analysis in which all 
stability roots were computed. It is emphasized that the approximate 
analysis is only approximate in that the calculation of the less signi- 
ficant stability roots is omitted. Clearly, by a judicious mix of 
complete and approximate analyses the total elapsed computer time 
required for the definition of an optimum engine design can be dras- 
tically reduced. 

4.1 THE FIRST TWENTY ITERATIONS 

Following initial experimentation using all stability roots , the engine 
was subjected to twenty design iterations using all specific combustor 
characteristics. An approximate analysis mode was employed which con- 
sidered only the relatively rapid calculation for the longitudinal time 
lag analysis (f 26 ) and the transverse time lag analysis (fpf) for the 
single Sv^ value of 3.0543 (the least stable transverse time lag root). 
This approximate analysis permits both longitudinal and transverse 
stability characteristics to be monitored. 

Initial and final stability roots from this optimization calculation 
are presented in Table VI. It can be seen that little change has 
occurred in the stability roots. The trend is to increased stability 
in the less stable transverse mode and to less stability in the more 
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stable longitudinal node. It may be noted that the transverse tine 

lag analysis of Table VI considers tvo solutions to the stability 

equation with S v = 3.05^3. These are the true solution with the 

freouencv approximating S,,. , and a spurious solution with the frequency 
- v n 

approximating r,/2. These spurious solutions with a frequency approx- 
imating tt/ 2 ;/e often encountered in the time lag analysis. If the 
true solution is not obtained on the nominal engine evaluation and the 
spurio'.s solution is obtained, the AUTQCOM program will "ti'ack" the 
spurious root. Hence, the analyst must take care to insure that the 
correct roots are found on the nominal design before embarking on an 
optimization run. This point is also discussed in Section 3.3. 


The engine rating after twenty design perturbations and the corres- 
ponding average engine characteristics are presented in Table VII. 

It can be seen that based on the selected average engine characteristic 
weights which provide the rating in the form of payload lost, a gain 
cf 25h pounds payload has resulted when compared to the nominal design 
cf Table III. It can also be seen that the average stability charac- 
teristic contribution to the rating is now negligible and that payload 
is being gained primarily by reduction of the performance characteristic 
penalty. Pursuing this payload improvement. Table VIII, it can be seen 
that the performance improvement stems from fpp » P e1 ' cen ^ mass of fuel 
vaporized, and from slight improvement in C* efficiencies for both pulsed 
and non-pulsed combustors. 

-.2 THE FIRST KUIJDR5D ITERATIONS 

Following the first twenty design iterations discussed in Section 4.1, 
the optimization problem was restarted without any stability analysis; 
and 100 successive design perturbations were introduced. A combination 
of the uniform random ray and pattern searches were employed. 
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F 1 1 

= 

100. 

F 1 2 

= 

100.000 

FI 3 


100.000 

fia 

= 

<51 .7069 

Fib 


73.6393 



— 

" 0 : “f 

F2 0 

= 

K21 

= 

0. 

F22 

= 

0. 

F23 

= 

0, 

F24 

= 

0. 

F2S 

= 

0. 

F26 

= 

-1638.16 

F27 

= 

-S68.3AS 

F2b 

= 

0. 

F29 

= 

°* - i 

F31 


87.9479 

F 32 

= 

50.0986 

F*1 

= 

433.943 

F42 

= 

10.0228 

F*3 

s 

9.97392 

F44 

= 

.399630 

FA5 

= 

6.005603E 

F<*6 

= 

Not Comp! 

FS1 

- 

2.20692 

F t> 1 

= 

2.08429 

F71 

s 

5 . 05579 


TABLE VIII. SPECIFIC CO’3USTOR CHARACTERISTICS AFTER 
TWENTY DESIGN PERTURBATIONS 
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Only f 2 g and f 2 j stability 
roots were computed 



I 


Reference 1. The approximate analysis employed completely neglects 
the stability characteristic. The rationale for this approach 
was the negligible stability characteristic contribution to the engine 
rating. Table VII. This table indicates that the stability charac- 
teristic affects the rating in the twenty- fourth significant figure. 

This is well below the accuracy of the CDC 6600 computer which, with 
sixty bits, is able to provide approximately ten significant decimal 
f i gure s . 

The nominal engine rating without the penalty of all stability charac- 
teristics (U.7 pounds. Table III) is 293.1 pounds. After 100 successive 
design perturbations introduced through the References 1 and 2 multi- 
variable search program, AESOP, the rating is reduced to 210.7 pounds. 
Rating convergence is illustrated in Figure 2. Convergence behavior 
of the combustor design variables is illustrated in Figures 3a through 
3c. The combustor design variables were allowed to fluctuate by plus 
or minus twenty- five per cent of the nominal values in this study. Two 
of the design variables, the chamber diameter and the number of fuel 
orifices (which equals the number of oxidizer orifices) are practically 
on the lower and upper bounds permitted in the study. 

The final rating and the characteristic components to the rating are 
presented in Table IX. Final design variable values together with 
the search limits employed are tabulated in Table X. From Tables III 
and IX the rating changes associated with each characteristic are 
seen to be 

Performance Characteristic 20.9 lbs., gain 

Stability Characteristic Wot considered 

Pressure Drop Characteristic O.lU lbs., gain 

Injector Complexity Characteristic 0.72 lbs., loss 


29 



FIGURE 2. COMBUSTOR RATING CONVERGENCE 
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TABLE X. DESIGN VARIABLE VALUES AFTER 100 PERTURBATIONS 



Length Characteristic 

Chamber Diameter Characteristic 

Chamber Mixture Ratio Characteristic 


36 . 7 lbs . , gain 

26 . 7 lbs . , gain 
1. 32 lbs . , loss 

Total Gain 82. U lbs. 

The total rating gain of 82. U lbs. produced by the optimization process 
of the AUTOCOM code ignores any stability characteristic effect. To 
assess this effect, a complete analysis was performed using the Table X 
vector of combustor design variables. The rating resulting from this 
complete analysis is presented in Table XI. The associated specific 
combustor characteristics are presented in Table XII. The stability 
characteristic produces a rating component of .16 pounds, a U.55 pound 
improvement over the nominal engine stability characteristic. Comparing 
the final rating of 210.8^ pounds. Table XI, with the complete nominal 
engine rating of 297-78 pounds. Table III, the total rating gain obtained 
in 100 design perturbations is 86 . 9 ^* pounds. It is interesting to note 
that despite the use of an approximate analysis which resulted in the 
stability characteristic being ignored, this characteristic nonetheless 
improved during the 100 design iterations. Elapsed computer time for 
the 100 iterations , the final complete analysis , and the initial complete 
analysis was 250 seconds on the CDC 6600 computer. 

It. 3 A NOTE ON STABILITY ROOTS AFTER 100 ITERATIONS 

The complete stability root set obtained after 100 iterations is presented 
in Table XIII. It can be seen that the second frequency corresponding 
to S v ^ - 3.05^3 is missing. This root was the least stable on the 
nominal engine. Table IV, but became more stable in the first 20 iter- 
ations of Section l+.l, Table VI. Accordingly, a search for this root 
was initiated to confirm the stability improvement over 100 iterations. 

The root was located as a non-conjugate pair at the points 
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TABLE XII. SPECIFIC CCMBUSTOR CHARACTERISTICS AFTER 
100 DESIGN PERTURBATIONS 
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TABLE XIII. STABILITY ROOT SET AFTER 100 ITERATIONS 


= (_.M* 3 09U + J2. 73775) 


and z 2 = (-.552853 + j3. 61289) 

Both roots are well damped; however, since the imaginary parts of these 
roots differs markedly from the value of 8^(3.051*3) a ’’ray search" was 
carried out through the design space. This search proceeded along the 
ray joining the nominal engine design to the final design obtained after 
100 iterations. The ability to carry out this type of ray search through 
an n-dimensional space (in this case, a twelve-dimensional space) is a 
standard feature of the AESOP program. Fifty-two points were equi- 
distributed along the ray search joining the nominal and final design. 

The root corresponding to S v ^ = 3.051*3 was tracked along the ray > starting 
from the nominal design. Root variation along the ray is presented in 
Figure 1*. The root at 

z = (-.03257 + J3.0521) 

presented in Table VI tracks continually into the root at 
z x = (-.1*1*3091* + j2. 73775) 

confirming this root as a valid solution to the stability root charac- 
teristic equation. Both final roots, zq and zp> obtained for S v ^ = 3.051*3 
are, therefore, considered to be valid roots. Their heavily dampeu 
nature results in their providing no contribution to the final engine 
rating. It can be seen from Figure It that the root at z^ is becoming 
more stable as the design progresses and that the root at Z2 is becoming 
less stable. 
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kA VERIFICATION OF THE OPTIMAL SOLUTION 

The optimal solution reported in Section 4.2 was verified in two ways. 
Firsts the solution was continued for 100 additional iterations with 
the uniform, random ray and pattern search algorithm. A slight perfor- 
mance improvement resulted. A final rating of 209.43 pounds was 
attained, a gain of 1.3 pounds over the solution of Section 2.2. 

Secondj the solution was restarted from the nominal solution using 
a different search algorithm. The algorithm used in this second 
solution was a recently developed directed random ray search, Appendix 
B, in c.mbination with the pattern acceleration algorithm. The final 
rating obtained by this method was 209.46 pounds after 100 iterations. 
Convergence of this solution has been added to Figure 2. It is clear 
that a final solution has been obtained. It is also clear that the 
newly developed search provides more rapid convergence to the solution 
than the older uniform directed ray search. This behavior is in keeping 
with other tests of the new search. 
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5.0 CONCLUSION 


The AUTOCOM code has successfully developed an improved engine design 
starting from the existing nominal engine. The payload potential of 
the engine was improved by 8? pounds as measured by the rating equation. 
Computer time required by the AUTOCOM code was minimal. The average 
time requirement for an assessment of each combustor design was approx- 
imately two seconds on the CDC 6600 computer. Computer time absorbed 
by the optimization subprogram AESOP in determining suitable design 
variable perturbations was negligible — approximately 103 seconds. The 
engine was optimized in one hundred design perturbations; hence, total 
computer time required to optimize the design was approximately four 
(1 ) minutes. More computer time would be required if combustor stability 
problems had been encountered. In this eventuality, it is estimated 
that twenty (20) minutes computer time would be required to obtain a 
solution. A definitive assessment of computer time in such a case 
awaits further experience using the AUTOCOM code. 

An examination of the optimal engine components reveals that the 
payload gain was largely obtained from improvements in the perfor- 
mance, chamber length, and chamber diameter characteristics. Small 
payload gains also resulted from improved stability and pressure 
drop characteristics. The injector complexity characteristic and the 
chamber mixture ratio characteristics both contributed performance 
losses when the final engine is compared to the nominal engine. 

A complicated set of design variable perturbations were introduced 
to obtain the payload capability improvement. An assessment of the 
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design variable changes introduced by the optimization algorithms 
indicates that the number of fuel and oxidizer holes, volume of the 
oxidizer dome, volume of the fuel dome, length of the combustion 
chamber, chamber diameter, and the mixture ratio are all sensitive 
design variables in the engine considered. In particular, in both 
optimal solutions obtained the number of fuel and oxidizer holes 
rapidly rises to the upper limit permitted, indicating that further 
payload improvement might result from a further increase in the 
number of holes allowed. Diameter of the fuel holes, diameter of 
the oxidizer holes, length of the fuel holes, and length of the 
oxidizer holes were relatively insensitive design variables for the 
engine design considered, presumably because of the basic stability 
of this engine. 



6.0 REFERENCES 


Hague, D. S. and Glatt, C. R.: An Introduction to Multivariable 
Search Techniques for Parameter Optimization (and Program AESOP). 
NASA CR-73200, April 1968. 

Hague, D. S. and Glatt, C. R. : A Guide to the Automated 

Engineering and Scientific Optimization Program - AESOP. NASA 
CR-73201, April 1968. 



APPENDIX A 

OUTLINE OF THE AUTOCOM PROGRAM ANALYSIS 


The AUTOCOM program automatically determines the combustor chamber 
characteristics given the chamber design variable values. The 
analysis considers performance, stability, and injector complexity 
characteristics. In an optional mode of operation, the program 
possesses the ability to automatically perturb the design parameters 
defining the engine characteristics (optimization). Stability and 
performance analysis modules available within the program are described 
below. 

Al, PER CENT MASS OF FUEL VAPORIZED (FUNCTION fn ) 

Per cent mass of fuel vaporized is computed by the method of NASA TR-67, 
Reference Al. 

"A model and theory for describing the rocket combustion 
process are described. The model is based on the 
assumption that propellant vaporization is the rate- 
controlling combustion process. Calculations of the 
vaporization rate and histories show the effects of 
propellants, spray conditions, engine design parameters 
and operating parameters on the vaporization process. 

The results are correlated with an effective chamber 
length for ease in using them for design purposes. An 
analysis is presented on the quantitative effect of 
incomplete propellant vaporization on combustor perfor- 
mance. With this analysis, experimental and calculated 
combustor performances are compared for injectors where 
drop size can be calculated. For other injectors the 
drop sizes are deduced and are shown as functions of 
injector type and orifice size." 


Al 



A2. PER CENT MSS OF OXIDIZER VAPORIZED (FUNCTION fi o) 

Per cent mass of oxidizer vaporized is computed by the method of 
Reference A1 in a similar manner to the fuel vaporization method 
summarized above. 

A3. C* EFFICIENCY BY MIXING MODEL (FUNCTION fi g) 

The first method available for computing C* efficiency is based on 
the method of NASA TN-2881, Reference A 2: 

"A model for predicting rocket combustion performance is 
presented which is based on the assumption that performance 
is limited only by gas-phase turbulent diffusion, or 
mixing, of oxidant and fuel vapors. The model shows how 
mixture ratio, chamber length, injector-hole spacing, 
and turbulence intensity affect performance. 

"Many physical processes occur simultaneously in a rocket 
combustor. In order to understand the importance of the 
various processes, such as vaporization, gas-phase mixing, 
or chemical reaction, each one is considered separately so 
that their effects on combustor performance may be deter- 
mined and compared. The vaporization process in rocket 
combustion is well understood, and an exhaustive analysis 
of it has been presented in the literature. Chemical 
reaction rates are usually considered to be very rapid and, 
therefore , not a limiting factor in controlling the rocket 
combustion process. A treatment of the relative importance 
of chemical reaction rates in rocket combustion is presented 
in Reference A3. The mixing process, though less understood, 
may possibly be, under certain conditions, a rate limiting 
step in the combustion process. 
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"In essence, the model developed in Reference A2 combines 
the highly generalized results of Reference Ah vith a 
technique similar to that suggested in Reference Al. In 
Reference Al, it is suggested that the effect of mixing 
on performance may be determined by calculating the perfor- 
mance of many small areas in a combustor cross section and 
averaging the results. The results of Reference Ah show 
how propellant concentration varies radially across the 
combustor as a function of chamber length, injector hole 
spacing, and intensity of turbulence, but do not indicate 
what effect such variations might have on combustor performance. 

"The method of NASA TN-2881 translates the generalized 
concentration profiles of Reference Ah into combustor 
performance values. A model based on that of Reference Ah 
enables mixing-limited performance to be calculated for 
particular propellant systems as a function of chamber length, 
turbulence intensity, injector-hole spacing, and operating 
propellant mixture ratio. Results of detailed digital computer 
calculations using this model are presented in Reference A2 
for eight propellant systems: oxygen with hydrogen, ammonia, 

hydrazine, and JP-h; fluorine with hydrogen, ammonia, and 
hydrazine; and nitrogen tetroxide with hydrazine." 

Ah. C* EFFICIENCY FOR PULSED COMBUSTORS (FUNCTION fi U) 

C* efficiency for pulsed combustors is computed by the statistical 
relationships presented in NASA CR-72370, Reference A5. 

"The objective of this method is the establishment of 
criteria for the design of stably operating liquid 
propellant rocket engines by means of a systematic 
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analysis of existing test data. In this analysis , 
relationships were sought between engine design 
variables, operating variables, and stability charac- 
teristics. The results of theoretical and experimental 
studies of combustion instability were used as guides 
in seeking these relationships . 

"The method was established by 

1. Development of a system for collecting 
rocket engine stability test data and 
utilization of this system to collect 
such data from a wide variety of engines • 

2. Definition and evaluation of functions 
of engine variables (parameters) which 

may be related to stability characteristics. 

3. Establishment of relationships between 
engine design and stability parameters by 
analysis of the collected experimental 
data. 

h. Formulation of an approach for utilizing 

these design-stability relationships in the 
development of new engines . 

"The results provide a comprehensive description of past 
experience with combustion instability in various engine 
types. The design approach offers a means for utilizing 
this experience to avoid development of new engines which 
are prone to instability." 



A5. C* EFFICIENCY OF NON-PULSED COMBUSTORS (FUNCTION fi g) 

C* efficiency for non-pulsed combustors is computed by the statistical 
relationships presented in NASA CR-72370, Reference A5 , discussed above. 

a6. chugging decay rate rased on fuel system (FUNCTION fo 0 ) 

Function f^ measures the fuel system chugging decay rate based on the 
method of Reference A 6 for either pump or pressure fed systems. Pump 
fed systems decay rates are found from the eigenvalues of the charac- 
teristic equation 2.04.03 of Reference A6. 

F(s ) = (1 + Es +JEs 2 ) [l + s - n + ne“ Ts ]+ PEse _TS = 0 (2.04.03) 

Pressure fed system decay rates are found from the eigenvalues of the 
characteristic equation (2.05.02) of Reference A6. 

F(s ) = [l + Js + JEys 2 + J 2 Ey (l-y )s 3 ] (l + s - n + ne“ TS ) 

+ Pe** TS (1 + JEys 2 ) = 0 (2.05.02) 

Eigenvalues are found by application of modern optimization procedures 
to minimization of |f(s)| followed by root sweeping. 

A7. CHUGGING DECAY RATE BASED ON OXIDIZER SYSTEM (FUNCTION f en) 

Function f 2 i measures the oxidizer system chugging decay rate by the 
method of Section A6 above. 

A8. STABILITY CHARACTERISTICS FOR PULSED OPERATION BASED ON STATISTICAL 
CORRELATION (FUNCTION To o) 

Function f ^ is the characteristic for pulsed operation based on the 
regression analysis of Reference A5. Basis of this approach is described 
in Section A4. 
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A9. STAB I LITY CHARACTERISTIC FOR NON-PULSED OPERATION BASED OW 
STATISTICAL CORRELATION (FUNCTION f 23 ) 

Function f^ is the characteristic for non-pulsed operation based on 
the regression analysis of Reference A5. Basis of the approach is 
described in Section Ah. 

A10. FULL SYSTEM HIGH FREQUENCY STABILITY DECAY RATE BASED ON THE 
MV ilQD OF DYKEMA (FUNCTION fp h) 

Function fph is the fuel system high frequency stability decay rate 
based on the method of Dykema, Reference A7- The characteristic 
decay rates of selected longitudinal and transverse mode combinations 
are computed. The Dykema method provides 

"A simplified engineering approach to the analysis of 
high frequency combustion instability in large liquid 
rocket engines. The approach stems from theoretical 
consideration of pressure and time dependent droplet 
combustion. There results a dimensionless correlating 
parameter called a stability number (N s ) which essen- 
tially represents the dimensionless ratio of a charac- 
teristic molecular diffusion time to a characteristic 
acoustic time. Stable and unstable ranges of N s are 
defined, and N s is reduced and simplified to common, 
readily measurable engineering terms involving the in- 
jector orifice pattern (size and number of orifices), 
the frequency of the acoustic modes, chamber pressure, 
and propellant flow rate . " 

All. OXIDIZER HIGH FREQUENCY STABILITY DECAY RATE BASED ON THE METHOD 
OF DYKEMA (FUNCTION f 

Function f ^ is the high frequency stability decay rate based on the 
method of Dykema, Reference A7. The characteristic decay rates of 
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selected combinations of longitudinal and transverse modes are 
computed. The Dykema method is summarized in Section A10. 


A12, STABILITY DECAY KATE PA SEP OR SENSITIVE TIME LAG MODEL FOR 
LONGITUDINAL MODE (TUKCTIK] f ?fi) 

Function f ^ is the stability decay rate characteristic based on 
the sensitive time lag model for a longitudinal mode. Reference A 6. 
The correlation equations for the; interaction index developed by 
Reardon of Aerojet is incorporated in the model. Decay rates are 
based on the eigenvalues of Equation (3-01.20) of Reference A6. 

— = M[(l - Yn) + yne~ TS ] (3.01.20) 

l+Be 25 

The solution is subdivided into 

a. Non-hyper golic propellant with coaxial injection 

b. Non-hype rgolic propellant with non-coaxial injection 

c. Storable propellants 


A13. STABILITY DECAY RATE EASED ON SENSITIVE TIME LAG MODEL FOR 
TRANSVERSE MODE (FUNCTION f 2 j± 

Function f ^ is the stability decay rate characteristic based on the 
sensitive time lag model of Reference A8. Decay rates are based on 
modifications to the characteristic Equations (28) of Reference A8. 


hjP + h 2 = 0 


\ = Yujl-js^ E f (u/u e )dz] 

Jo 
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— ytN 

s . nz. 



h 0 = -(Y + 1)5. - j(f - 7)s v>1 Ze 

z e 

j (u/u g )dz] 

0 

a. Non-hype rgolic propellant vith coaxial injection 

b. Non-hype rgolic propellant with non-coaxial injection 

c. Storable propellants 

All*. STABILITY DECAY RATE BASED ON THE RESPONSE FUNCTION APPROACH OF 
LEWIS RESEARCH CENTER (FUNCTION f ? ft) 

Function is the decay rate determined from the acoustic wave solutions 
of Priem and Rice, Reference A9- Response functions for liquid propellants 
are determined by Reference A10. Response functions for gaseous propellants 
are determined by Reference All. 

A15. STABILITY CHARACTERISTIC BASED ON THE NON-LINEAR ANALYSIS OF 
PRIEM_ AND GUENTERT (FUNCTION f oo) 

Function f ' is the decay rate based on the non-linear analysis of Priem 

and Guentert , Reference A12. 

"Regions of combustion instability in rockets are calculated 
from a non-linear theory that considered the combustor to 
be an annular section with very small thickness and length. 

Two models are used to determine the local burning rate. 

One assumes that the burning rate is equal to the vapori- 
zation rate; the other assumes that the burning rate is 
equal to the chemical-reaction rate. The results show that 
a finite disturbance is required to produce instability. 

The instability regions are found to be a function of 
several design parameters and to be insensitive to the 


C ± LL 

+E £ (*•-?> + ^ + i)s„ h 5 e 


As in Function f 26 the model is specialized for 
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activation energy, specific-heat ratio, and order of 
reaction of the propellants. The vaporization rate 
motel is more sensitive to a pressure disturbance for 
design parameters corresponding to conditions encountered 
in large combustors. The chemical-reaction-rate model 
is more sensitive to a pressure disturbance for conditions 
corresponding to small research combustors. Wave shapes 
and characteristics are determined for various conditions. 11 

Al6. engine design and complexity characteristics 

The remaining functions in Figure A1 are straightforward engine design 
and complexity factors. 


f^ is the fuel pressure drop characteristic 


f 2 is the oxidizer pressure drop characteristic 


f^ is the number of fuel plus oxidizer holes characteristic 


f^ is the volume of the oxidizer dome characteristic 


f^ is the volume of the fuel dome characteristic 


fj^ is the length of the oxidizer holes characteristic 


is the length of the fuel holes characteristic 


f » r is an injector-type complexity characteristic 


f is the chamber length characteristic 
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fgi is the chamber di are ter characteristic 


is the mixture ratio characteristic 
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APPENDIX B 

DIRECTED RANDOM RAY SEARCH 


This search proceeds along a succession of random rays distributed 
about a best estimate of the gradient vector. The search can be 
used in combination with the pattern search acceleration procedure 
of the AESOP program. 

The best gradient vector estimate, R, is based on a weighted combi- 
nation of the old gradient vector estimate, *^e latest 

search step direction which improved performance, R f . 

R i = (W R ’ Ri oia + R i' )/(W R + 1 *°) (El) 

The search step directions explored are based on a weighted combi- 
nation of the best gradient vector estimate, R, and a small random 
vector, r. 


6a. = (U R • B i + r.)/(U R +'l) (B2) 

On problems involving a pronounced ridge in the control space, this 
search will prove efficient. Once the approximate direction of the 
ridge is established by a performance improvement, the random rays 
are focused in the general direction of the ridge, and excursions 
outside the region of improvement tend to be minimized. The search 
is sensitive to the weighting constant values Up and Up. Based on 
a study of the Rozenbrock Valley problem, nominal values of Wr = 5«0 
and Ur = 2.5 are recommended. 
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It should be noted that while the directed random ray search proves 
effective when an approximate ridge direction is defined. Figure Bl, 
it may prove wholly ineffective when the ridge abruptly changes 
direction, Figure B2, or when acquisition or a ridge requires a 
large directional change , Figure B3. To avoid convergence failure 
in these last two situations, the weighting constant, U r , which 
focuses the random rays must be adaptively determined. When further 
progress proves impossible for a given value of Up, this weighting 
constant must be decreased. As U r -*■ 0, the search approaches the 
uniform random ray search which permits an abrupt change of search 
direction. Following establishment of a new search direction, the 
random rays are refocused along the new approximate ridge direction 
by an increase in Up. Logic to focus and defocus the directed 
random rays is included in the AESOP code. 
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APPENDIX C 


WEIGHTING FACTOR CONSTANTS USED FOR 
THE SAMPLE CASE OF THE 15,000 LBF ENGINE 


Ayj = 150.0 Constant of the performance characteristic in 

the rating equation. 

Apjj = 200.0 Constant of the stability characteristic in 

the rating equation. 

^FIII = .00069 Const ant of the pressure drop characteristic 

in the rating equation. 

Apiv = -0395 Constant of the injector complexity characteristic 

in the rating equation. 

Ayy = U 0.5 Constant of the chamber length characteristic 

Apyi = 35*0 Constant of the chamber diameter characteristic 

Ayvn = 66.0 Constant of the mixture ratio characteristic in 

1 ~ the rating equation. 


B = 1.0 Exponential on the performance characteristic 

in the rating equation. 

Bp^ = 0.0 Exponential on the stability characteristic in 

the rating equation. 

Bfih = 1.8 Exponential on the pressure drop characteristic 

in the rating equation. 

Bp^y =2.0 Exponential on the injector complexity charac- 

1 teristic in the rating equation. 

= 1.86 Exponential on the chamber length characteristic 

Bpyj =2.0 Exponential on the chamber diameter characteristic 


Cl 


FVII 
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Exponential on the mixture ratio characteristic 
in' the rating equation. 



Exponential on the stability characteristic 
in the rating equation. 


fll 


0.0 


Constant in the performance characteristic 
equation. 


a fl2 


1.0 


Constant in the performance characteristic 
equation. 


fl3 


1.0 


Constant in the performance characteristic 
equation. 


a fii+ 


.01 


Constant in the performance characteristic 
equation . 


fl5 


.01 


Constant in the performance characteristic 
equation. 


af 20 = 1.0 


a f21 ' 1 '° 


a f22 " 1,0 


a f23 * °-° 


a f2lj " 0 '° 


a f25 1 * ° 


8 f26 * 1 -° 


a „_= 1 .0 

f27 


Constant in the combustor stability characteristic 
eqxmtion. 

Constant in the combustor stability characteristic 
equation. 

Constant in the combustor stability characteristic 
equation. 

Constant in the combustor stability characteristic 
equation . 

Constant in the combustor stability characteristic 
equation. 

Constant in the combustor stability characteristic 
equation. 

Constant in the combustor stability characteristic 
equation . 

Constant in the combustor stability characteristic 
equation. 
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Constant in the combustor stability characteristic 
equation. 

Constant in the combustor stability characteristic 
equation . 

Constant in the fuel pressure drop characteristic 
equation. 

Constant in the oxidizer pressure drop characteristic 
equation 

Constant for the injector orifice number. 

Constant for the oxidizer dome volume. 

Constant for the fuel dome volume 

Constant for the length of the oxidizer orifices. 

Constant for the length of the fuel orifices. 

Constant for the injector type complexity. 

Constant for the chamber length characteristic 
Exponential on the chamber length characteristic 

Constant for the chamber diameter characteristic 
Exponential on the chamber diameter characteristic 

Constant for the mixture ratio characteristic. 

Constant for the mixture ratio characteristic. 

Exponential on the mixture ratio characteristic 
equation. 

C3 



